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Abstract. We consider the basic model of virus dynamics in the modeling of 
Human Immunodeficiency Virus (HIV), in a 2D heterogenous environment. It 
consists of two ODEs for the non-infected and infected CD^ T-lymphocytes, T 
and /, and a parabolic PDE for the virus V . We define a new parameter Aq as 
an eigenvalue of some Sturm-Liouville problem, which takes the heterogenous 
reproductive ratio into account. For Ao < the trivial non-infected solution is 
the only equilibrium. When Aq > 0, the former becomes unstable whereas there 
is only one positive infected equilibrium. Considering the model as a dynamical 
system, we prove the existence of a universal attractor. Finally, in the case of 
an alternating structure of viral sources, we define a homogenized limiting 
environment. The latter justifies the classical approach via ODE systems. 

1. Introduction. The acute infection by the Human Immunodeficiency Virus (in 
short HIV) is characterized by a huge depletion of the CD^ T-lymphocytes {CD4) 
and a peak of the virus load [7]. After few weeks, these two components reach a 
steady state which characterizes the asymptomatic phase of the infection. Before 
the availability of highly active antiretroviral therapy, this later phase lasted after 10 
years in median with an accelerated decrease of CD4 and an increase of virus load. 
A substantial number of nonlinear ODE systems have been suggested by Perelson 
et al. (see [21, 22]) to model the complex dynamics of HIV-host interaction. For 
instance, such models have been used to estimate the infected cell half-life and 
the viral clearance during antiretroviral therapy [11, 20, 29], or to understand the 
dynamics during acute infection [23]. 
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The common basic model of viral dynamics [2] includes three variables: T, the 
non-infected CD4, I, the infected CD4, and V, the free virus: 

Tt ^ a ^ jVT ~ tiTT, (1.1a) 

It^"fVT^^lII, (1.1b) 

Vt = Nfiil - fivV. (1.1c) 

This model, describing the interaction between the replicating virus of HIV and host 
cells, is based on some simple hypotheses. Non-infected target cells are produced by 
the thymus at a constant rate a and die at a rate /xt- By contact with the free virus 
particles (virions) they become infected at a rate proportional to their abundance, 
'yVT (see Fig. 1). These infected cells die at a rate m and produce free viruses 
during their life-time at a rate N. Free particles are removed at a rate called the 
clearance nv- AH these parameters are generally positive constants. This simple 
model study has led to interesting results (see [18, 28]) and suggested a treatment 
strategy (see [2]). 



2- HIV uses CD^ 
machinery to make 
more virus 




Figure 1. The HIV life cycle 



It is easily seen that the system has two equilibria: 

(i) the non-infected steady state 

Tu = — , lu^ 0, K = 0, 

/iT 

which corresponds to a non-negative equilibrium in case of no infection; 

(ii) the infected steady state 

7A' 7A' 7/iy 

also called seropositivity steady state, corresponding to a positive equilibrium 
in case of infection. 

Some authors (see e.g., [2, 18]) have considered the basic reproductive ratio Rq: 

Ro = (1.3) 
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a dimensionless parameter defined by epidemiologists as the average number of 
infected cells that derive from any one infected cell in the beginning of the infection 
[18, p. 16]. Stability properties of the two steady states are usually studied around 
this quantity: if i?o < 1 the non-infected steady state is stable, if i?o > 1 the 
infected steady state has a biological meaning and it is stable, while at _Ro = 1 
both steady states coincide. So, i?o = 1 is a bifurcation point (see Fig. 2). Further 




Figure 2. Typical bifurcation diagram, with Rq varying between 
and 3, as a function of a. The solid line represents the stable branch 
of V, the dashed line the unstable one. The other parameters have 
the values: 7 = 0.001, N ^ 1000, fir = 0.1, /i/ = 0.5, ^lv = 10. 
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Figure 3. Profiles of virus (left) and target cells (right) in the 
case of infection. Here i?o = 1-5 and the other parameters have the 
same values as in Fig. 2. 

models have been used involving other populations present in the immune system 
(see [17, 18, 4]). However, these models assume that the populations T, I, V are 
homogeneous over the space for all time, which is a common, but not a very realistic, 
assumption. Actually, the interaction between the virus and the immune system 
(either as a target with CD4 or as an agent for controlling infection) is localized 
according to the type of tissues [3] and also in a given tissue (e.g. lymph nodes). 
To examine the effects of both diffusion and spatial heterogeneity, Funk et al. [8] 
introduced a discrete model based on (l.la)-(l.lc). These authors adopted a two- 
dimensional square grid with 21 x 21 sites and assumed that the virus can move 
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to the eight nearest neighboring sites. They pointed out that the presence of a 
spatial structure enhances population stability with respect to non-spatial models. 
However, our analysis does not confirm this observation (see Section 4 below). 

Recently, Wang et al. [26] generalized Funk et al.'s model. They assumed that 
the hcpatocytcs can not move under normal conditions and neglected their mobility, 
while viruses can move freely and their motion follows a Fickian diffusion. They 
proposed the following system of two ODEs coupled with a parabolic PDE for the 
virus: 

Tt^a- jVT - ^itT, (1.4a) 

It = iVT - fill, (1.4b) 

Vt^ NfijI ^ fivV + dvAV, (1.4c) 

where dy is the diffusion coefficient. They assumed that the domain is the whole 
real line and proved the existence of traveling waves. Wang et al. [27] introduced 
a delay to take into account the time between infection of a target cell and the 
emission of viral particles [6]. They considered (1.4a)-(1.4c) in a one-dimensional 
interval with Neumann boundary conditions. 

In the spirit of the above works, we intend to study System (1.4a)-(1.4c) in a 
two-dimensional spatial domain (0,^) x (0,£) with periodic boundary conditions. 
There are two main situations: 

(i) the environment is homogeneous and, hence, all the parameters in (1.4a)-(1.4c) 
are constant. Therefore, the system with diffusion has the same equilibria as 
System (l.la)-(l.lc); 

(ii) the environment is heterogeneous, therefore certain parameters become posi- 
tive functions of the space variable. Then, the virus is spatially structured. 

For simplicity, we assume throughout the paper that only the rate a varies while 
the other parameters arc fixed positive constants. In fact, it is biologically plausible 
to assume that the arrival of new CD^ may vary according to local areas. More 
precisely, a is piecewise continuous and periodic in each variable with period £. 
Then it is convenient to define the heterogeneous reproductive ratio: 

i?o(x) = 2iX^, x e (0, £) X (0, i). (1.5) 

The sites where Ro{x) < 1 are called sinks while the sites where Ro{x) > 1 are 
called sources [8]. 

The paper is organized as follows. In Section 2 we are interested in the stationary 
problem associated with (1.4a)-(1.4c) and its non-negative equilibria. The virus 
equilibrium verifies the elliptic semilinear equation 

V 

dv^V - HvV ^ -fJ.Tt^vRo{x)—ry-^ , (1.6) 

with periodic boundary conditions. A first issue is to define a parameter which will 
play the role of the bifurcation parameter i?o in the case the latter is constant. A 
candidate for this role is the largest eigenvalue Ao of the operator: 

dvl^ + fiv{Ra-l)Id, (1.7) 

which is the linearization around V = oi (1.6). For the reader's convenience, we 
recall some basic facts about two-dimensional Sturm-Liouville eigenvalue problems 
with periodic boundary conditions such as (1.7) and give some proofs in Appendix 
A. 
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It is clear that, whenever Rq is a constant, Ao = fiv{Ro ~ !)• Therefore, we 
distinguish two cases, depending upon the sign of Aq: 

(i) -^0 < 0: the trivial non-infected solution = is the only solution of (1.6); 

(ii) Aq > 0: (1.6) has exactly two non-negative solutions, namely the trivial non- 
infected solution Vu = and the positive infected solution Vi. 

Section 3 is devoted to the study of the evolution problem (1.4a)-(1.4c). In the 
case Ao < 0, we prove that the trivial non-infected solution {Ti,Ii, Vi) is asymptot- 
ically stable. Then, we turn our attention to the biologically relevant case Aq > 0. 
First, we prove the non-infected solution becomes unstable. Second, we consider 
(1.4a)-(1.4c) as a dynamical system and prove the existence of an universal (or 
maximal) attractor. Since the system is only partly dissipative, we use a result of 
Marion [16]. The following Section 4 is devoted to some special cases where the 
positive infected solution Vi is stable. Particular attention is paid to the case when 
Rq is a constant: in this case discrete Fourier transform can be applied. 

We point out that our proof can be extended to further models in HIV litera- 
ture. It is not difficult to take a logistic term into account in the T equation [22], 
although the steady equation (1.6) will be more involved. Adding such a term, Hopf 
bifurcations have been observed numerically in ODE systems (see [20]). Therefore, 
proving the stability of the infected solution may be, in general, challenging. 

In the last section (Section 5), we consider the case when a heterogeneous envi- 
ronment is formed of sinks and sources alternating very rapidly, with a heteroge- 
neous reproductive ratio i?o(f). We determine the homogenized limiting medium as 
e — > 0. It is fully characterized by a constant reproductive ratio, the mean value of 
Rq. Therefore, the classical approach of HIV dynamics via ODEs in a homogenous 
environment can be a posteriori justified in this respect. 

Notation. Throughout this paper, for any > 0, we denote by the usual space 
of functions / : (0,^)^ M such that is integrable. The square (0,£)^ will be 
simply denoted by fJ^. By H'' we denote the Sobolev space of order k, i.e., the 
subset of of all the functions whose distributional derivatives up to k-th order 
are in . Both and are endowed with their Euclidean norm. Finally, by 
we denote the closure in H'^ of the space C™ of all m-th continuously differentiable 
functions / : ^ K which are periodic with period £ in each variable. The space 
is endowed with the norm of H'^. Finally, we denote by Id the identity operator. 

2. A semilinear equation for the virus steady states: existence and unique- 
ness of the equilibria. Wc start from the system for the virus dynamics: 



set in fl£. Periodic boundary conditions for T, / and V are prescribed. 

We are interested in the existence of steady state solutions to the equations 
(2.1a)-(2.1c) which belong to the space x x H^. Clearly, any steady state 
solution to Problem (2.1a)-(2.1c) is a solution to the following stationary system: 



Ii 



a — "fVT — ^itT, 

^VT - fiji, 

NfijI ~ i^vV + dvAV, 



(2.1a) 
(2.1b) 
(2.1c) 



a — 'yVT — ijtT = 0, 
-/VT - inl = 0, 



(2.2a) 
(2.2b) 
(2.2c) 
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From a biological point of view, only non-negative solutions to (2.2a)-(2.2c) have a 
meaning. Hence, we limit ourselves to proving the existence of this kind of steady 
state solutions. 

System (2.2a)-(2.2c) can be reduced to a single scalar equation for the unknown 
V. Actually, it is not difficult to infer from (2.2a), (2.2b) that 



7y + ^T' /i/(7y + ^T)' 

Hence, the function V turns out to solve the equation 

^aNV V 
dvAl/-^yl/ = ^ -^lT^lvRo-r^, , (2.3) 

associated with periodic boundary conditions. 

2.1. Existence and uniqueness of non- negative equilibria. In this subsection 
we will provide a thorough study of the equation (2.3). As it has been already 
stressed, we are interested in non-negative solutions only. 

Clearly, equation (2.3) always admits the trivial non-infected solution = and, 
hence. Problem (2.1a)-(2.1c) admits 

T„(x) = ^, /„(x)=0, K(a;)=0, a; e 17,, (2.4) 

as a (trivial) steady state solution. We will call the triplet {TutIu.Vu) the non- 
infected solution. 

We are interested in studying the uniqueness of the non-infected solution in the 
class of all the non-negative steady state solutions to Problem (2.1a)-(2.1c). Of 
course, in the case when uniqueness does not hold (a situation which can actually 
occur, look for instance at the case when i?o > 1 and a is constant, discussed in the 
introduction) we want to characterize all biological relevant steady state solutions 
to Problem (2.1a)-(2.1c). 

For this purpose, we need to recall the following results about Sturm-Liouville 
eigenvalue problems in dimension two with periodic boundary conditions. 

Theorem 2.1. Let d and /i be, respectively, a positive constant and a bounded 
measurable function. Further, let £/ : Lr' be the operator defined by £/u — 

dAu — fj.u for any u G . Then, the spectrum of si consists of eigenvalues only. 
Moreover, its maximum eigenvalue Amax is given by the following formula: 

A„.ax = - inf I ' • (2.5) 

Finally, the eigenspace corresponding to the eigenvalue Amax is one dimensional and 
contains functions which do not change sign inQi. 

This is a rather classical result. Nevertheless, for the reader's convenience, we 
give a proof in Appendix A. 

In view of Theorem 2.1, we can define the constant Aq to be the maximum 
eigenvalue of the operator ip i-^ dyAtp -\- fiv{Ra ~ l)<y5, which is the linearization 
around V = of operator V y-^ dyAV — ^vV + ^Tl^vRoiy^^:^- According to (2.5), 

. r f dv Ja^iipx^dx + fiv JnS^ - Ro)4'^dx\ 
- Ao = mf <^ jj—^ ^ . (2.6) 
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As we are going to show, the uniqueness of the non-infected steady state solution 
is related to the value of Aq. 

Lemma 2.2. // Ao < 0, then the non-infected solution {Tu-, Iu,Vu) {see (2.4)) is 
the only non-negative solution of (2.1a)-(2.1c). 

Proof. We argue by contradiction. Let us suppose that Problem (2.2a)-(2.2c) admits 
another solution {T*,I*,V*) different from {Tu,Iu,Vu)- Then, the function V* € 
does not identically vanish in and it solves the equation (2.3). Multiplying 
both the sides of this equation by V* and integrating by parts in ilf, we get: 

dy f \VVrdx + t,vf )ivrdx = 0, 

J a, JnA lV*+iiTj 

or, equivalently, 

dv I |Vy*pdj; + ^y I {I- Ro){V*fdx 

+ ^iv I RoU- ) {Vn'dx = 0. 

Jn, \ iV +1^tJ 

Since V* does not identically vanish in Qi, the last integral term is positive, implying 
that 

dv f \VV*\'^dx + fiv [ (1- Ro){V*)^dx <0. 

Hence, the infimum in (2.5) is negative which contradicts our assumption — Aq > 
0. □ 

Remark 2.3. From formula (2.C) it is immediate to check that, when the maximum 
of i?o in is less than or equal to 1, the constant Aq is non-positive. Hence, in 
this situation the non-infected solution is the only relevant steady state solution 
to Problem (2.1a)-(2.1c) in complete agreement with the case when Rq is constant 
(see the Introduction). 

The result in Lemma 2.2 is very sharp as the following theorem shows. 

Theorem 2.4. Suppose that Aq > 0. Then, there exists a steady state solution 
{Ti,Ii,Vi) to Problem (2.1a)-(2.1c) whose components are all positive in fif, with 

= ^ , I, = 1^ . (2.7) 

Moreover, {Tu,Iu, Kt) 0,'rid {Ti,Ii, Vi) are the only steady state solutions whose com- 
ponents are non-negative in fl/. 

Proof. It is clear, that we can limit ourselves to dealing with the equation (2.3). 
Being rather long, we split the proof into two steps. 

Step 1: {existence). To prove the existence of a positive solution to the equation 
(2.3) in H"^, we use the classical method of upper and lower solutions. To simplify 
the notation, we denote by ^ the sup- norm of the function Rq. We look for an 
upper solution vq of (2.3) as a constant C > 0. It is immediate to check that the 
best choice of wq is 

vo{x) = , X £ 

7 

Note that, by Remark 2.3, M is strictly greater than 1. 
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To determine a lower solution, in the spirit of [14, Chapt. 13, Sec. 3], wc take as a 
candidate to be a lower solution the function Vq = cipQ with c > to be fixed. Here, 
tpo is the unique (positive) solution to the equation dyA(/3o + l^v{Ro ^ ^)fo ^ ^o'i^o 
which satisfies sup^^gj^ ^0(2;) = 1- 

If we plug cipo in (2.3), we get 

dvA{c(po) - nvcfo + fJ-TfJ-vRo ^ 

= XoC(po + fiTf^vRo — fivRoClfiO 

japQ + fiT 

= cipo [ Ao - 7Mv^o — 

V -yctpo + ^iT 

Mr / ■ 

Since Ao > is fixed, the last side of the previous chain of inequalities is non-negative 
as soon as Ao/iT ~ ^c^v^ ^ 0- Hence, if we fix 

. r Ao mt(^-I)' 
c = mm < , 

the function ~ cipo turns out to be a positive lower solution to the equation (2.3) 
and it satisfies < vq in fie. 

Hence, the classical method of upper and lower solutions provides us with a 
positive solution to the equation (2.3). It is enough to define the sequence (u„) by 
recurrence in following way: for any fixed n G N, TJ„ is the unique solution in of 
the equation 

dvAVn - fiVVn + HT fJ-V Ro _ ^ = 0. 

7W„_1 + /.JT 

Since the function 1 1— > h{t) := is increasing in [0, +00), the maximum princi- 

ple (see Proposition B.l) shows that the sequence (jJn) is pointwise non- increasing. 
Moreover, Uq <Vn < vq for any n G N. Hence, the sequence (U„) is bounded. More- 
over, by very general results for the heat equation, there exist positive constants Ci 
and C2, independent of n, such that 

WVuWh^ < Ci||i?o^(l^n-l)||L2 < C2||iJ„_i||l2, 

tt 

for any n G N. Thus, the sequence (w„) is bounded in as well. Since it converges 
pointwise in we can now infer that converges strongly in and weakly in 
to a function v G which, of course, turns out to be a solution to the equation 
(2.3). For further details on the method of lower and upper solutions, we refer the 
reader, e.g., to the monograph [19]. 

Step 2: (uniqueness). To prove the uniqueness of the nontrivial non- negative 
solution to the equation (2.3), we adapt to our situation a method due to H.B. 
Keller [13]. 

Let us suppose that u is another nontrivial non-negative solution to (2.3). Then, 
the function w := v ^ u belongs to and solves the equation 

dyAw ~ fiyw + fi^fivRoqiu,v)w = 0, (2.8) 

where 



Qix, y) = ; r- ■ -, x,y>Q. 

(7.T + fj.T)[iy + fJ-T) 
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Let US denote by Arnax(u,w) and Amax(w, 0), the maximum eigenvalues in of the 
operators 

di/A + ^^^vRoq{u, v)Id 

and 

dy A + fi^^vRoq{u, 0)ld, 
respectively. By Theorem 2.1, they are given by the formula 

. , [dyL iVV'l^dx-^l/iy L i?og(w,jw)'0^rfa;l 

for j = 0, 1. Since v is non-negative and it does not identically vanish, q(u,v) < 
g(it, 0) and there exists an open subset of Qi where q{u,v) < q{u,0). Hence, 
Amax(u,w) < Amax('", 0). Clcarly, since w satisfies (2.8) and it docs not identically 
vanish in Qe, then < Amax(M, v) < Aniax('u, 0). 

Let us now rewrite the equation satisfied by u in the following way: 

dyAu - Aniax(M, 0)u + ^|^yi?o9(u, 0)u = (^y - Ainax(w, 0))u If. 

Fredholm alternative implies that ip should be orthogonal to C, where by C we have 
denoted the function which spans the eigenspace associated with the eigenvalue 
AmaxCujO). As it has been already remarked (see Theorem 2.1), the function ( 
does not change sign in fie. Similarly, since nv < Knax{u, 0), (/s is non-positive and 
it does not identically vanish since u docs not. Hence, the function (p cannot be 
orthogonal to C and this leads us to a contradiction. The proof is now complete. □ 

2.2. Numerical illustration (steady state). In accordance with Funk et al. [8], 
the domain Cli is a discrete square grid with nxn sites of equal dimension £/n x £/n. 
We assume in this numerical part that all parameters vary randomly from site to 
site in such a way that 0.1 < Ro{x) < 5.0. We deal with two cases: 
(i) In the first case, see Fig. 5 (left), the distribution of Rq is as in Fig. 4a. The 
sources represent only 26% of the sites. Wc compute Aq = —1.70. Solving 
numerically the equation (2.3), we find the non- infected solution Vu = 0. We 
also represent T„ according to Formula (2.4). 




Figure 4. Two distributions of Rq on a 40 x 40 grid, (a): 26% of 
the sites are sources; (b): 50% of the sites are sources. 




Figure 5. Densities of virus V (top) and target cells T (bottom). 
Left: Rq as in Fig. 4a, Aq < (no infection). Right: Rq as in Fig. 
4b, Ao > (infection). Here dy = 1, £ = 1, N = 40. 



(ii) In the second case, see Fig. 5 (right), the distribution of Rq is as in Fig. 
4b. Now the sources represent half of the sites. The eigenvalue Aq = 4.30 is 
positive. Numerically, we observe the positive infected solution Vi of (2.3). 
Note that Vi is smoothly structured in space although Rq is not. We also 
represent Ti according to Formula (2.7). 

3. Study of the dynamical system. We recall the evolution problem for the 
virus dynamics: 

Tt^a--fVT-fiTT, (3.1a) 
It='fVT-niI, (3.1b) 
Vt^ Nfiil ~ HvV + dvAV, (3.1c) 

set in rii with periodic boundary conditions. We consider (3.1a)-(3.1c) as a dy- 
namical system ^{t), which has two equilibria: the non- infected trivial solution 
(Tti, Iu,Vu) and the infected, positive solution (T^, Ii,Vi), the latter for Aq > only. 
At first, we prove that the non- infected solution is stable for Aq < and unstable 
for Aq > 0. By stable we mean asymptotically stable. For Aq > the instability of 
(Tk, Iu,Vu) does not usually imply the stability of (Ti, Ii,Vi). Our aim is to prove 
the existence of a universal (or maximal) attractor which attracts all the orbits (see 
e.g., [25]). Since System (3.1a)-(3.1c) is only partly dissipative, we will use a result 
of Marion [16]. Some special cases where the stability of the infected solution is 
granted will be discussed afterwards. 
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Let US introduce the following notations: '3 is the domain in K"^ defined by: 

•2) = {{x, y,z) eM.^ -.0 <x + y < Mi, 0<z< A/2}, 

where Mi, M2 arc positive constants which will be fixed throughout the proof of 
the next theorem. We also set for u := (wi, U2, U3): 

H = {u G (L2)3 : u{x) e 3 for a.e. x e Qi}, V = x x . 

Our main result is the following theorem. 

Theorem 3.1. The following properties are met: 

(i) if Xq < 0, the trivial non-infected solution to Problem (3.1a)-(3.1c) is stable; 

(ii) if Xq > 0, the trivial non-infected solution to Problem (3.1a)-(3.1c) is unstable; 

(iii) if Xq > and a € W^'°°{fli), the dynamical system S^{t) associated with 
(3. la) -(3.1c) possesses a universal attractor that is connected in H. 

3.1. Proof of (i). We begin the proof observing that the linearization (around 
(T„, /„, Vu)) of Problem (3.1a)-(3.1c) is associated with the linear operator ^„ de- 
fined by 



^ —pLxId 


V 





-fijld 



70; 
7a 



Id \ 



Id 



Nfiild dyl^-pLvIdJ 

Its realization L„ in [L^Y with domain D{Lu) = x x generates an ana- 
lytic strongly continuous semigroup. Indeed, is a bounded perturbation of the 
diagonal operator 





/-pirld 








\ dvA-fivIdJ 

defined in x x H^, which is clearly sectorial since all its entries are. Hence, we 
can apply [15, Prop. 2.4.1(i)] and conclude that L„ is sectorial. Since is dense 
in L^, the associated analytic semigroup is strongly continuous. 

Let us prove that all the elements of the spectrum of Lu have negative real part. 
In view of the linearized stability principle (see e.g., [10, Chapt. 5, Cor. 5.1.6]) this 
will imply that the trivial non-infected solution to Problem (3.1a)-(3.1c) is stable. 

To study the spectrum of the operator Lu, we fix f = ifi, f2, fa) S (i^)^ and 
consider the resolvent system 

(A + ^t)^i + — '/'s^/i, (3.2a) 
(A + M/)<P2-— V3 = /2, (3.2b) 

{X + fiv)V3 - dv^V3 ~ N^ii(p2 = fs, (3.2c) 

where we look for a triplet of functions <pi,</?2 G and ip^ S H^. Suppose that 
A differs from both —^j and —fxr (which belong to the essential spectrum). Then, 
we can use equations (3.2a) and (3.2b) to make tpi and ip2 explicit in terms of 
ip3. In particular, replacing the expression of ip2 in terms of ip^ into (3.2c) and 
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using the very definition of the function Rq (see (1-5)), we can transform Problem 
(3.2a)-(3.2c) into the equivalent equation for (p^ only: 

(A + iiv)'P3 ~ dvAtps - !^\^^ Roipz = --— + h- (3.3) 
A + /i/ A + /i/ 

Adding and subtracting ~iivRq from the left-hand side of (3.3), we can rewrite the 
equation (3.3) into the equivalent form: 

(A + ^y)<^3 - dyA(/33 - fivRom + ^ivRoT-^ — V3 = + h- (3-4) 

A + /i/ A + /i/ 

Note that, for any A G C, the operator A\ defined by the left-hand side of (3.4) 
has compact resolvent. Hence, its spectrum consists of eigenvalues only. We are 
going to prove that, for A G C with non-negative real part, is not an eigenvalue 
of Ax- For this purpose, we observe that — Aq can be equivalently characterized as 
the infimum of the ratio 

dv In, \'^i^\^dx + fiv Ja^{l - RoM'^dx 

when tp runs in the set of all the complex- valued functions tjj = tjji + iijj2, with 
V'iiV'2 G Hj. This shows, in particular, that 

-Xof \tp\^dx<dv f \\7ip\^dx + fiv f {l-RoM\^dx, (3.5) 

for any function ?/; as above. 

Let now (^3 be a complex- valued solution to (3.4). Multiplying both the sides of 
such an equation by the conjugate of (f^ and integrating by parts, we easily see that 



(X + fiv) / l^psl dx + dv I iVi^sl dx - ^.v / i?o|</'3| dx 
Jf, 

My A 



4- 



/ Ro\ip3\''dx = 0. (3.6) 
Taking the real part of both the sides of (3.6) and using (3.5), we obtain 
(ReA + ^y) / W-i\'dx-\J \^^\^ dx + ^Jiv^^^^^^ I W^?dx <Q. 

JUi Jill l^ + M-fl Jcte 

Since, by assumptions, Aq < and i?o is a positive- valued function, the only solution 
to the previous inequality, when ReA > 0, is the trivial function (^3 = 0. Hence, 
for these values of A, is not an eigenvalue of the operator A\, i.e., any A with 
non- negative real part belongs to the resolvent set of the operator L„. 

3.2. Proof of (ii). Again in view of the linearized stability principle, to prove the 
instability of non- infected solution (T„,/„,l^) we can limit ourselves to showing 
that Lu admits a positive eigenvalue. Hence, we are led to the study of Problem 
(3.2a)-(3.2c) with /i = /2 = /s = 0. Since the parameters /it, and /iy are all 
positive and we are looking for positive eigenvalues A, we can limit ourselves, as in 
the proof of (i), to studying the equation 

dyA(/33 + I /f£^ _ tA ^ 

VA + /i/ J 
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For any s > 0, let us consider the operator dyA + iiy — 1^ Id defined in 

H^. By Theorem 2.1 its spectrum consists of eigenvalues only, and the largest one 
is given by the following formula: 

dy J^^ \Vifdx + Hv (l - f^) ^^dx ■ 



^A(,s) = inf 



s+fJ.1 



As it is immediately seen, A(0) = Aq. Moreover, since the function s 
is continuous, decreasing in [0, +oo) and it tends to as s ^ +oo, the function 
s I— > A(s) is continuous, decreasing and tends to 

A(+cxj) = - mf < — } , 

as s ^ +00, i.e., it converges to the largest eigenvalue of the operator dy A — ^y/d, 
which clearly is — /xy. Now, since s i— > A(s) is a decreasing continuous function 
mapping [0, +(X)) into [— /.ty, /Xy^Ao], it is immediate to check that the fixed point 
equation s ~ A(s) has a positive solution. Of course, this fixed point is the positive 
eigenvalue of the operator L„ we were looking for. This accomplishes the proof. 

3.3. Proof of (iii). Let us show that [16, Thm. 5.1] applies. In this respect the 
assumption a e W^'°°{^li) is enough. To avoid conflict with notations, throughout 
the proof, T denotes time as usual, whereas the triplet (T, /, V) is denoted by 
M2, U3). We split the proof into several steps. 
Step 1. Here, we prove that, for any Uq := (uq,!, uo,2j 1*0,3) G i^'^Y 1 the Cauchy 
problem 

Dtui{t, •) = -HTUi{t, ■) - jui{t, ■)u3{t, ■) + «(•), t > 0, 
DtU2{t, •) = -mu2{t, ■) + 7Ui(t, ■)\uz{t, •)!, ^ > 0, 

(3.7) 

Aw3(t, •) = Au3(i, •) - A*yW3(i, •) + NniU2{t, •), t > 0, 

Ui(0, ■) = Mo,i, j = 1, 2, 3, 

Problem (3.7) admits a unique classical solution defined in some time domain 
(0,r*). Here, by classical solution, we mean a vector valued function u such that 
ui,U2 e Ci((0,r.);L2) nC([0,T,);L2) and U3 G C([0, T.); L^) n C7i ((0, T,); L^) n 
C((0,T,);i/2). 

Problem (3.7) is semilinear with a nonlinear term which is a continuous function 
from X L2 X DA2{f3, 2) into L"^ x L"^ x L"^ for any (i G (1/2, 1). Here, D/^^[I3, 2) is 
the interpolation space of order (/3, 2) between and the domain of the realization 
A2 of the Laplacian with periodic boundary conditions in D{/S.2) = H^). 

Hence, D/^^{I3,2) = {L'^,H^)f3^2 and this latter space coincides with H^^ , which 
continuously embeds into the space Q of all continuous and periodic (with period 
i in each variable) functions (see e.g., [9, Thm. 1.4.4.1]). 

To prove the existence of a classical solution to Problem (3.7), let us fix /3 G 
(1/2, 1) and introduce, for any T > 0, the space Cfj{T) consisting of all functions 
V : (0,r) Da2{P,oo) such that ||t;||c^(T) := suptg(o_-r] |li^w(i, ■)\\D^^if3,oo) < +00. 
Clearly, Cp{T) is a Banach space when endowed with the above norm. Moreover, 
D\^{f3, Qo) is continuously embedded into D\^{P — e,2) for any e G (0,/?). We 
now fix e small enough such that 6 = (3 — e > 1/2. From the above results, it is 
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immediate to infer that Ci3{T) is embedded into the set of all continuous functions 
/ : (0, T] X — > R and there exists a positive constant Ci, independent of T, such 
that 

sup t''||/(t,-)l|oo<Ci||/||c,(T). (3.8) 

te(o,T] 

Let us solve the Cauchy problem for ui, taking as a parameter. The (unique) 
solution to such a problem in ((0, T]; i^) nC([0, T]; L^) is the function ui defined 

by 

iti(i, •) =exp ^-^Ti - 7 j U3{s,-)ds^ uo,i 

+ a{-) J exp ^-/^T(t - s) - 7^ U3{r,-)dr^ ds, (3.9) 

for any t G (0,T]. If uo,i is bounded and continuous in Qi this result is straight- 
forward. In the general case, we approximate uo.i G by a sequence of smooth 
functions Uq'i. It is immediate to check that the function it^"' defined by (3.9), 
with Uq"]' instead of uo.i, converges to the function ui in (3.9) in C([0,T];L^), by 

dominated convergence. Similarly, Dtu"i^ converges to Dtui in C([e, T — e]; i^) 
for any e > 0. It follows that the function ui in (3.9) is a solution to the Cauchy 
problem for ui also in the case when tio,i is in L^. 

Let us now denote by Ai the operator defined in Ci3{T) by the right-hand side of 
(3.9). A very easy computation shows that Ai maps Cp{T) into C^"^([0, T]; X^). 
Moreover, 

llAi(«)||c([o,T];L^) <exp(^^7Ti"'ll«llc,(T)) (ho,i||L^+T||a|U2), (3.10a) 
||Ai(u) - Ai(u;)||c([o,T];L2) 

+ (l_g^(2-g) ^^'''"^'i^-) II"' - "llc^MT)- (3.10b) 

We now consider the equation for U2. Replacing ui — Ai(w3) in the right-hand 
side of this equation and using the same argument as above, we easily see that the 
(unique) solution in C([0, T]; L^) n ((0, T]; L^) the function U2 defined by 

U2{t,-) = e-*'"uo,2 + l I e-'^^(*-^V3(s,-)l(Ai(w3))(s,-)ds, [0,T]. (3.11) 



Let us denote by A2 the operator defined in Ci3{T) by the right-hand side of 
(3.11). Taking (3.8), (3.10a) and (3.10b) into account, one can easily show that 

C 

l|A2(w)||c([0,T];L2) < ||mo,2||l2 + J^"''^^^'' \Mc,3{T) (II"o,i||l2 + T||a||L2) 

xexp('^7Ti-''Mlc,(T)), (3.12a) 



\\My)\\c^-Oi[0,T]-L-) < (^^' + 1) II"0^2||l = 



HETEROGENEOUS VIRAL ENVIRONMENT IN A HIV SPATIAL MODEL 15 

+ Y^7(l + T'-')\\v\\c,iT) exp (^^jT'-'^MciT) 

xi\\uo,i\\L-+T\\a\\L2), (3.12b) 

||A2(u;) - A2(u)||c([0,T]:L2) 

X jhoallL^ + tPWl^ + \\v\\c,iT) (i + (T^^^^^'^'ll^^ll^^) } ' 

(3.12c) 

for any v,w G Cf}{T). Let us now observe that (ui,M2,U3) is a classical solution to 
Problem (3.7) if and only if U3 is a fixed point of the operator F, formally defined 

by 

(r(«))(i,.) = e-*'^^e*^=?/o,3+ / e-^-(*-^'e(*-^'^HA2(i;))(s,-)rfs, tG(0,r). 



We are going to prove that the operator F is a contraction in ^^(r) = {it G 
Cp{T) : ||u||c3(T) ^ provided that T and M are properly chosen. As a first 
step, let us prove that F maps ^p{T) into itself if T,M are suitably chosen. For 
this purpose, we set 

A-=suppV^^|b,^(^,^). 
t>o 

Taking (3.12a) into account, we can estimate 

l|r(z;)||c,(T) < jlh^o.slU^ + Y^ho,2||L= 

I rT.2-e _____ f C"! 

for any v G '3^p{T). Hence, if we fix M > K\\uo,3\\l'^, we can then choose T small 
enough such that ||F(i')||c^(j') < M for any v G ?Vp{T). Moreover, taking (3.12c) 
into account, we can estimate 



(Tr^^'exp (^^^ri-«Af) (||«o,illc,(T) +T||a|U.)| 



||F(H -r(tO||c,:,(T) < ^T||A2(U;) - A2(«)|ic([0,T];L^) 

< ^r^^f^'-' exp II. - .11,,., 

X |ho,illL^ + T^II^IU^ +m(i + J^^T^J^^'T'^'Ml^ 

for any VjW G '3^g{T). This estimate shows that F is a 1/2-contraction provided 
that T is sufficiently small. We can thus apply the Banach fixed point theorem 
and conclude that there exist T > and a unique function G '3'p(T) solving the 
equation F(u3) = u^. 

The function us actually belongs to C([0, T]; L^) n C\{0, T]; L"^) D C((0, T];Hf). 
Indeed, by (3.12b), the function A2(u3) is in C^-^{[0,T]; L'^). Therefore, [15, Thm. 
4.3. l(i)] guarantees that the function F(m3) has the claimed regularity properties. 
Moreover, Dtu^ = Au3 — /iyU3 + A2(it3) in (0,T]. As a byproduct, the triplet 
(iti, 1*2,^3) is a classical solution to Problem (3.7). 
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By a classical argument we can extend the solution (ui,U2,M3) to a maximal 
solution defined in some time domain [0, T*). This vector valued function (still 
denoted by (lii, ^2, U3)) enjoys the following properties: wi,U2 S C([0, r,);L^) fl 
Ci((0, T^);L^), U3 G C([0, T,); L^) n ^^((0, T,)-L^) n C((0, T,); H^). 

Step 2. Here, we prove that if Uq > (where the inequality is meant componen- 
twise) then the maximal defined solution u to Problem (3.7) is non- negative as well 
in (0,T*). Clearly, using formulae (3.9) and (3.11) it is immediate to check that ui 
and U2 are both non-negative whenever uq^i is. 

Let us now consider the problem for W3, which we rewrite here: 

r A"3(t,-) = Au3(t,-)-My"3(i,-) + ^M/«2(i,-), iG(0,r,), 

{ (3.13) 

[ "3(0, ■) = Uo,3- 

The heat semigroup is positive in Cj by the maximum principle. Since the heat 
semigroup in Cj is the restriction to Cj of the heat semigroup {e*^^} in , by 
density it follows that {e*'^^} is non-negative as well. This is enough for our aims. 
Indeed, the function U3 is given by the variation of constants formula 



W3(t,-)=e*('^^-^-)7/o,3 + iVM/ / e(*-^)('^^-^-)u2(s,-)ds, tG(0,T,), (3.14) 



uo,3 and U2 being non-negative, the function U3 is non-negative as well. 

We have so proved that any solution to Problem (3.7), corresponding to an initial 
datum in the first octant, is confined to the first octant for any t £ (0,T,t). In such 
a case we can forget the absolute value in (3.7). 

Step 3. Here, we prove that any solution u to Problem (3.7), corresponding to 
an initial datum Ug £ exists for any positive time and it stays bounded. Here, 

Si = {(x,y,z) eM.^ ■.Q<x + y< A/i, < z < M2}, 

ll«ll°o j,,T MiNfii 

A/i = 7, M2 = ll"lloo- (3.15) 

mmj^T,/^/} t^v 

For this purpose, it is convenient to introduce the so-called Svab-Zeldovich variable 

V ~ ui + U2. As it is immediately seen, the function v satisfies the Cauchy problem 

f Dtv(t,-) = a - fiTUi - fiiU2, te(0,T*), 

I v{0, •) = Uo,l + "0,2- 

Since ui and U2 are both positive, then 

Dtv{t, •) < a -min{^T,M/}'"(^> ■) -.^ a - ^iv{t, ■) < ||q;||oo - t^-vit, ■), t £ (0,T*). 

(3.16) 

Multiplying both the sides of (3.16) by a non-negative function (p G L'^ and integrat- 
ing over one obtains that the function w{t) = J^^ v{t, x)ip{x)dx is in C^([0, T,)) 
and solves the differential inequality 

Dtw{t) < \\a\\^\\ip\\Li(^Q^) ~ fiw{t), te (0,T,). 

Hence, 



w{t) < e-''* f / v{0,x)ip{x)dx + i(e^* - l)||a||oo / ip{x)dx) , t G (0,Th.), 
or, equivalently, 

/ (v{t,x)~e-^''v{0,x)--il-e-^'*')\\a\\^)ip{x)dx<0, t e {0,n). 
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From this integral inequality, we can infer that 

v{t,x) < e-^'v{0,^) + -i'^ - e~^')\\a\\oo < -\\a\\^ ^ Ah, (3.17) 
M A* 

for any t G [0, T,) and almost any x G il^. 

Since v ~ ui + U2 and ui and U2 arc both non-ncgativc, it follows that ui and 
U2 can be estimated by the right-hand side of (3.17). 

Finally, let us consider the function u^. From (3.14) and the above results, we 
can infer that 

DtU3{t,-) < Ausit,-) ^ fivuait,-) + MiNfii, t e (0,T,). 

Note that the function 1*3 = A/2, with M2 being given by (3.15), satisfies the 
previous inequality. Hence, if 7/9,3 M2, then, the solution to Problem (3.13) is 
bounded from above by M2. With the previous choices of Mi and M2, we see that 
the solution to Problem (3.7) which corresponds to uq G stays in & for any 
t G (0, r*). By virtue of [15, Prop. 7.1.8], u can be extended to all the positive 
times. 

Step 4- Here, we show that, for any Uq G the solution u that we have 
determined in the previous steps is, in fact, the unique weak solution to Problem 
(3.7) which belongs to L'^{{0,T)- L^) x L^{{0,T)- L^) x L^{{0,T); Hi) for any T > 0. 
Even if the following arguments are standard, for the reader's convenience we go 
into details. 

As a first step, we observe that, since mi, U2 and U3 are bounded, the weak 
derivatives _Df Ml and DtM2 are in i°°((0, +00) x fif ) = L°°((0, 00); i°°(r2f )). Hence, 
ui and U2 are locally Lipschitz continuous in [0, +00) with values in L°^(Qe). 

Let us now consider the Cauchy problem for U3 (i.e., problem (3.13)). Since U2 is 
Lipschitz continuous in [0, T] with values in L^, for any T > 0, by [15, Thm. 4.3.1(i)], 
such a Cauchy problem admits a solution which is in C([0, T] ; L^) n ((0, T] ; L^) n 
C((0,T]; H^). By the weak maximum principle, the Cauchy problem (3.13) admits 
a unique weak solution. Hence, U3 G C{[0,T];L^) n C^{{0,T];L^) n C{{0,T]; Hf). 
Now, we turn back to the equations for ui and U2 and conclude that DfUi and DtU2 
are in C([0, T];L^) for any T > 0, this implying that ui and U2 are in C^([0, T]; L^). 
Hence, any weak solution to Problem (3.7) with data in & is such that ui,U2 G 
C1([0,Th.);L2) and 1*3 G C{[0,T,); L^) D C\{0,T,); L^) H C{{0,n); Hf). Since we 
have proved uniqueness of the solution in this class of functions, uniqueness of the 
weak solution follows as well. 

Since for non-negative solutions the Cauchy problem (3.7) coincides with problem 
(3.1a)-(3.1c), we have, thus, established the following: 

Proposition 3.2. For every Uq G H, the Cauchy problem (3.7) possesses a unique 
solution for all time, u(t) G H for all t, u G L'^{0,T;V) for all T > 0. The 
mapping S^{t) := Uq i-^ u(t) is continuous in H. Furthermore, if Uo G V, then 

Step 3. We are now in a position to apply [16, Thm. 5.1]. The set 

& = {{x, y,z) eR^ -.0 <x + y < Mi, 0< z< M2}, 

is a positively convex, compact region of M"^. To meet all the hypotheses of [16, Thm. 
5.1], it remains to consider the non-dissipative part of (3.7), i.e. the equations for 
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Ml and U2 which we rewrite in the compact form: 

;,)(:)+(») ft)"'"- 

Obviously the matrix G(-,U3) has positive eigenvalues whenever U3 > 0, which are 
bounded from below by positive constants. Hence, 

for any ^ £ K^. Thus, condition (4.6) in [16] is satisfied. The proof of Theorem 3.1 
is completed. 

4. A gamut of some special cases. 

4.1. Numerical illustration (evolution). We continue the discussion of Subsec- 
tion 2.2 in the framework of the evolution problem (3.1a)-(3.1c). In the first case 
(Aq < 0), only the non-infected steady state Vu = exists. We solve (3.1a)-(3.1c) 
under particular initial conditions: we start the infection at the center of the grid 
with an inoculum of one viral unit, assuming that T and / are at their uninfected 
steady state. One observes that the virus vanishes very rapidly and the target 
cells return to their initial level (see Fig. 6 left) in accordance with the stability 
of the uninfected equilibrium Vu- In the second case (Aq > 0), two equilibria exist, 
Vu and Vi. Starting with the same initial conditions, the virus population grows 
while the population of target cells decreases. Both of them achieve an equilibrium 
corresponding to the positive infected solution (see Fig. 6 right). 




Figure 6. Dynamics of free virus (top) and target cells (bottom) 
on a 40 X 40 grid. Left: Rq as in Fig. 4a, Aq < (no infection); 
right: Rq as in Fig. 4b, Aq > (infection). The solid, the long and 
the short dashed lines depict respectively the dynamics at (10,10), 
at the site of viral inoculum (20,20) and at a border site (1,1). Here 
dv = 1, i=l. 
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We are now going to review some particular cases where the stabiUty of the 
infected solution is granted. 



4.2. Homogeneous environment with diffusion. We consider the case when a 
is a constant such that i?o (sec (1.3)) is a constant, as well, greater than 1, together 
with c?y > as in [8]. 

The linearization around {Ti,Ii,Vi) (see (1.2)) of Problem (3.1a)-(3.1c) is asso- 
ciated with the linear operator 



^ —^rRoId 







V 







The same arguments as in the proof of Theorem 3.1(i) show that the realization 
Li of the operator in [Li^Y with domain D{L.i) ^ x x generates an 
analytic strongly continuous semigroup. We are going to determine its spectrum. 
For this purpose wc use the discrete Fourier transform. 

Wc consider the realization A2 of the operator A with domain H^. Its real 
eigenvalues can be labeled as a non-increasing sequence (— Afe), fc = 0, 1, . . . Only 
Aq = is simple, the other eigenvalues being negative such that — — > —00 as 
k — > +00. 

We claim that the spectrum of the operator Li is given by 



fceNu{o} 



(4.1) 



where <7k is the spectrum of the matrix 

/ "MTi?0 
- ^t{Ro - 1) -^ij 

V N^lJ 



Mv 

N 

N 

-Xkdv — 



fiV/ 



To check the claim, let us observe that, if a function v = {vi,V2,V3) in L^xL^xH^ 
solves the resolvent equation Au — L^u = f , for some A G C and f = (/i, /2, /a) in 
(L^)^, then its Fourier coefficients (wi,fc, W2,fc, W3,fc) {k = 0, 1, . . .) solve the infinitely 
many equations 



( A + ^(Ti?o 

^j.t{Ro-'^) A + /.t/ 
\ ~Nni X + Xkdv+fivJ 



' N 



V2,k 



(fiA 

f2.k 

\h.k) 



0,1, 



where fi^k denotes the k-ih Fourier coefficient of the function fi {i ~ 1, 2, 3). Clearly, 
any eigenvalue of Mj. (fc = 0, 1, . . .) is an eigenvalue of Li. Therefore, 

fceNu{o} 
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On the other hand, if A ^ crfe for any k = 0,1,..., then all the coefficients 
(■ui.fc, ?^2,fc, I's.fc) are uniquely determined through the formulae 

vi.jc = ^^(^x) { ^ + '^'^^^ + ~ fJ'Wv]fi,k - /i/Mv/2,fc 

- ^(A + Ai/)/3,fc}, 

V2,k = ^^^^-^ |mt(-Ro - 1)(A + Afcrfy + Hv)fl,k + (A + /lT-Ro)(A + Afcdy + ^v)f2,k 
+ ^(MT + A)/3,,}, 

«3,fc {Nfli^iTiRo - l)/l,fe + AfAi/(A + MTi?o)/2,fc 



where 



and 



+ (A + AJTi?o)(A + Ai/)/3,fe}, 
^fc(A) = A3 + 6fcA' + CfeA + dfe=0, fc = 0,l. 



Ck = ^J■TRo{^^I + A*v) + dvXk{^J■TRo + m/), 

4 = M/MvMt(^o - 1) + dvXkfJ.ifJ.TRo- 

Note that, if A differs from both — /i/ and —fixRo, then 

^fc(A) ~ {A^ + \{fJTRo + fJi) + i-JifJTRo} Afe, as k 

Hence, for any A ^ {-m, -/-ItRq} U UfceNu{o} ^k, it holds that 

^ A + MJ ^ 

X'^ + X{fiTRa + fJi) + fJifJxRo ' 

Ht{R.Q ~ + (A + fJTR.Q)f2,k 



V2,k 



V3,k 



X^ + A(/iTi?o + fJi) + hiHtRo 

NlIifiT{Ro - l)fl,k + NfJiiX + flTRo)f2M + (A + fJTRo){X + fll)h,k 1 



dv{A2 + A(/.jTi?o + /i/) + M/MT-Ro} Afc' 

as fc +00. It follows that the sequences {f 1,^-}, {w2.fc} and {AfcVs^^} are in i"^ . This 
shows that the series whose Fourier coefficients are ui^fc, V2.k and wa^^, respectively, 
converge in (the first two ones) and in H"^ (the latter one). The inclusion 

(j{Li) C {-Aij, -Mt-Ro} U (J 0-/C 

*;eNu{0} 

follows. We now observe that a straightforward computation shows that A = —fixRo 
is in the essential spectrum of Li. Also A ~ — ^/ belongs to the essential spectrum 
of Li and, in the case when /i/ = fix, it belongs also to the point spectrum. The 
set equality (4.1) is proved. 

Clearly ak has three eigenvalues (counted with their multiplicity), either all real, 
or one real and two complex conjugates. Routh-Hurwitz criterion enables us to 
determine whether the elements of Uk have negative real parts. The latter holds if 
and only if bk, dk and bkCk — dk are positive, which is clearly true whenever i?o > 1. 
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Remark 4.1. As it is easily seen cro is the spectrum of the operator 
/ -firRoId -^Id\ 

\ Nfijld -fiyldj 

which is associated with the hnearization at {Ti,Ii,Vi) of Problem (l.la)-(l.lc). 
Since, as we have already remarked, the spectrum of Li is the union of the sets ak 
(fc G N U {0}) and the points — /x/, — /xt-Ro, the scenario is one of the following: 

(a) the diffusion does not improve the stability of the solution {Ti,Ii,Vi). Therefore, 
the stability issue is identical to that of the ODE system {dy ~ 0); 

(b) the diffusion worsen the stability of the solution (T,, Jj, V^). 

Therefore, we are unable to confirm Funk et al. [8], who pointed out that the 
presence of a spatial structure enhances population stability with respect to non- 
spatial models. Only some smoothing effect can be credited to the diffusion. 

4.3. Death rates fij < fiTt > 0. This case leads to a mathematically interesting 
framework although it has little biological relevance, since fij > in the literature 
(e.g., = 0.01, fij = 0.39 in [5]). For the latter reason we will not elaborate the 
case extensively. 

Using the Svab-Zeldovich variable i; = ui +M2 we can transform Problem (3.1a)- 
(3.1c) into the following equivalent one for the unknowns ui, v and u^: 

' Dtui{t,-) = -fiTUiitr) ~ iui{t,-)u3(t,-) + a{-), i > 0, 

Dtv{t, •) = a - fiiv{t, •) - {fiT - fii)ui{t, •), t > 0, 

< DtU3{t,-)^Au3{t,-)^fivU3it,-) + Nniv{t,-)~NfiiUi{t,-), t > Q, 

Ui{0, ■) = uo,i, « = 1,3, 

. u(0, •) = Uo.l + Uo-2. 

(4.2) 

It is not difficult to see that the mapping ui is non-increasing, so is the 

mapping ui i— > w thanks to the hypothesis jiT — ^ii > Q. Hence, the mapping ^ v 
is non-decreasing. Finally, the mapping M3 i— > M2 = w — ui is non-decreasing. Based 
on these observations, following [19] it is possible to construct two sets of monotone 
sequences which converge to the solution of (4.2). These sequences start respectively 
from upper and lower solutions defined as in Section 2, to stay away from the trivial 
solution. It is well-known that a solution of an evolution problem constructed via 
such a monotone sequence scheme, with suitable initial conditions between upper 
and lower solutions, achieves a stable equilibrium (see [19]). Therefore, the infected 
solution is asymptotically stable. Numerical computations in the phase plan (see 
Fig. 7) illustrate the difference in the virus dynamics when /i/ < fix (monotonicity) 
and III > fiT (spirals). 

4.4. Quasi-steady problem. In this part we assume that T and / arc at their 
equilibrium. In such a case, (3.1a)-(3.1c) reads 

= a — jVT — iitT, 

= - fill, 

Vt^ Nuil - fivV + dvAV, 
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Figure 7. Virus vs. target cells in the phase plane, in the case 
of infection. Left: mean values. Right: the solid, the long and 
the short dashed lines depict, respectively, the values at (10,10), at 
the site of viral inoculum (20,20), and at a border site (1,1). The 
parameters have the following values: 7 = 0.001, N = 1000, fir = 
0.1, nv = 10, dy = 1, i = 1; a varies between 0.1 and 5.0; /ij = 1 
in (a), m = 0.01 in (b). 



which is equivalent to the scalar parabolic equation for V only, with periodic bound- 
ary conditions: 

Vt = dvAV-fivV + fiTf^vRo-r7^ • (4.3) 

The latter is the natural evolution problem associated with (2.3). It is clear that 
(4.3) has the same non-positive equilibria, namely Vu ~ and Vi > 0, in the case 
when Ao > 0. The stability of Vi can be proved according to [10, Sec. 5.3] by 
constructing a Lyapunov function. 

5. Homogenization. This section is concerned with the case wherein the environ- 
ment is heterogeneous and is formed of rapidly alternating sinks and sources. For a 
fixed integer k, we imagine that il^ is divided into a network of periodic squares 
Q,^, where e — £/k. The heterogenous reproductive ratio Rq will depend upon e, 
see Fig. 8b. The idea of homogenization is to let e ^ and find the equivalent 
homogenized medium. Therefore, such a heterogenous environment can be replaced 
by its homogenized limit for easier computations and analysis. 
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Figure 8. Periodic structure with £ = 1 {a): fti is divided in a 
network of 40 x 40 sites. The sites are pieced together in 4 x 4 
periodic squares Qe with s = 0.1. (b): zoom of fJ^ which contains 
16 sites. The darkness represents the scale of Rq on each site (the 
darker the square, the larger is Rq). 



More precisely, we introduce a normalized periodic function Rq as a function of 
the variable y = (2/1,2/2), of period (1, 1), and we define: 

Rq{xi,X2) ^ Rq ( — , — 
\ e £ 

Remark 5.1. x = {xi,X2) is the macroscopic variable while y = (2/1,1/2) is the 
microscopic one. 

We consider the problem: 

dvAV ~ fivV = -^.TPivRo^T^ , (5.1) 

on Qi with periodic boundary conditions as above. The idea is to find the limiting 
homogenized equation as e ^ 0. Wc start with the following lemma: 

Lemma 5.2. Let Rq : — > R 6e a periodic [with period 1 in each variable), 
piecewise continuous function. For x G R^ and e > 0, set Rq{x) = i?o(f )• Then, 
the following properties are met. 

(i) Rq tends to ./#j,(i?o) as e — > 0, weakly in ^^^^(R^) for any p > 1; 

(ii) For any £ > 0, set 

[ ^v/ojVV'P^^ + Mv/o,(l- W^^ l ,,,, 
— Aq — mt < -5— > . lo-^J 

i,eHli.^o[ Jn.'^^dx J 

Then, Ag iiv{-^y{Ra) — 1) as £ — > 0. 
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Proof. Property (i) follows straightforwardly from e.g., [12, p. 5]. Thanks to (i), 
one can take the limit as e ^ in the right-hand side of (5.2) and show that Ag 
tends to the largest eigenvalue of 

dyA-^y(l-.^y(i?o))M 

on fl( with periodic boundary conditions. Let us prove this claim. As a first step, 
we observe that there exist two constants Ci and C2 such that Ci < Ag < C2, since 
the function Rq is bounded. 

Next, Aq is the largest eigenvalue of the Sturm-Liouville eigenvalue problem 

dvAip - ^y(l - i?g)V' = A?A, 

with periodic boundary conditions we documented in Theorem 2.1, associated with 
the cigcnfunction -0^. We may assume that J^^^{ip^)'^dx — 1. As we pointed it out 
in Theorem 2.1, ip^ does not change sign. 

It is clear that -0^ is bounded in H^. Then, there exists an infinitesimal sequence 
{e„} such that Ag" — ^ Aq, 0^" — + ip'^ weakly in H^, strongly in H'^^^ and (hence) 
uniformly in fi^. Note that J^-^^{ip'^)'^dx = 1 and does not change sign. 

Since 

it is not to difficult to pass to the limit in the above equation as n ^ +00 in the 
distributional sense and see that 

dyAV'°-My(l--^y(^o))^° = A!]V°, (5.3) 

with periodic boundary conditions. Therefore, Ag is an eigenvalue of the operator 
dyA — nv{^ — Ro)Id, associated with the eigenfunction ip'^. Since ip'^ does not 
change sign, Ag is the largest eigenvalue of (5.3). Obviously, tp'^ is a constant, hence 
Ag is explicit. Finally, checking that all the sequence Ag converges to Aq, as £ ^ 0, 
is an easy task. □ 



Next, we prove the following result. 

Theorem 5.3. Assume that ^y{RQ) > 1. Then, there exists Eq such that, for any 
e G (O,eo]j the equation (5.1) has a positive solution G H'^. As e ^G, tends 
to 

V° ^^i^yiRo)-!), 

1 

in and uniformly in ili. 

Proof. As a first step, we observe that, since ./#j,(i?o) > 1, from Lemma 5.2(ii) it 
follows that there exist 5, Eq > such that Xq > 6 for e £ (0,eo], and (5.1) has a 
unique positive solution V'^. Clearly, AV'^ is bounded in L°°{ni) and this implies 
that is bounded in H^. Arguing as in the proof of Lemma 5.2, one can extract 
an infinitesimal sequence {£«} such that V^" converges strongly in H^~^ to a limit 
G H"^, which verifies the equation 

dyAy° - /iyF° = -^,^^,yj^ [R^)^l . (5.4) 

7 1/" + 
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Because of the periodic boundary conditions, Equation (5.4) has only constant 
solutions. Therefore, 

-fiyV = -fJ,Tf^V-^y{Ro)- 



7^° + Mt' 

and it only remains to prove that V'^ is not the trivial solution. With obvious 
notations, we recall (see the proof of Theorem 2.4) that V'^ > Vq = c^f% where (y9g 
is the positive eigenfunction associated with the largest eigenvalue Ap and 

' — ^ > mm i 



IS — mill \ -— , ( ^ mill s 

Since max (ySg = 1, it is clear that remains bounded away from whenever e < £o- 
Finally, checking that itself converges to as e ^ is immediate. This 
concludes the proof. □ 

5.1. Numerical illustration (homogenization). We consider a model where i?o 
is as in Fig. 8, taking its values on the elementary 4x4 grid (b) as in Tab. 1. 



i 


i?o(l,j) 


i?o(2,j) 


i?o(3,j) 


i?o(4,j) 


1 


1.60 


1.41 


1.55 


0.819 


2 


0.800 


0.165 


2.59 


0.872 


3 


1.20 


0.489 


1.37 


0.453 


4 


2.09 


4.25E-4 


0.270 


2.80 



Table 1. Values of i?o(j,j), i^j = 1,...,4, on fig as in Fig. 8b, to 
be viewed as a 4 x 4 matrix. 



It is easy to compute the mean value ^y{Ro) = 1.16 and the homogenized viral 
density V'^ = 16.1. Fig. 9 shows how the virus at e = 0.1 (left) oscillates slightly 
around its homogenized limit (right). 




Figure 9. (left) vs. V° (right). The grid and parameters are 
as in Fig. 8 and Tab. 1; e = 0.1, dy = 1, i ^ 1. 
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Appendix A. Proof of Theorem 2.1. 

Proof. It is well known that the realization A2 of the Laplacian in , with domain 
H^, is a sectorial operator. Since is a bounded perturbation of the Laplacian, 
s!/ is sectorial as well, and, hence, its spectrum is not empty. Let us fix Aq G 
p{^) n (/Zoo, +00) such that [Ao,+oo) is in the resolvent set of T. Here, by /ioo 
we denote the sup-norm of the function yu. The operator ^ :~ — XqI turns 
out to be invertible and its resolvent set contains [0, +00). Since the operator T^^ 
is continuous from into H^, it is continuous, in particular, from Hj into itself, 
when this latter space is endowed with the inner product 

{v,w) ~ d / {DxivDxiW + Dx2vDx2w)dx + / {fi + Xo)vwdx 
:= d Vv ■ Vwdx + (m + Xct)v'wdx, 

J J 

which is equivalent to the Euclidean inner product of H^. Since is compactly 
embedded into Hj (see e.g., [1, Thm. 3.7]), the operator ^q"^ compact from 
into itself Moreover, is self-adjoint in . Indeed, 

{£/q'^u,v) =d / V£/(^^u ■ Vvdx + / + Xo)£/Q'^uvdx 



d I '^uvdx+ I (/i + Ao)^ ^uvdx 

uvdx, (A-.l) 



for any u,v & L^. Now, from the general theory of self-adjoint compact operators, 
it follows that the spectrum of consists of a sequence of real eigenvalues which 
converges to 0. As a byproduct, the spectrum of j^/q consists of a sequence of 
eigenvalues diverging to —00. More precisely, (t(^o) C (— cx),0) and A G cr(^o) if 
and only if A^^ is in (t(^^^). In particular, the maximum eigenvalue of £/o is the 
inverse of the minimum eigenvalue of Since is a compact operator, its 

minimum eigenvalue is defined by 



Amin(^o ^) = iiif 



Taking (A.l) into account, we can estimate 



Amin('2^o ^) = ipf 



In, ^^dx 



ipeH^ ^#0 I dv \yxip\'^dx + J^^in + Xo)ilj'^dx 



^(zH^^^^o \ dv Jn, I V^-iAPda; -I- {p. + X„)'ijj^dx J 

Formula (2.5) follows at once, observing that Amax(-2^) = Ao + Amaxl-sz^)- 

The last assertion of the theorem follows from the Krein-Rutman Theorem ap- 
plied to the restriction of jz/q to the space Ctt(R^) (of all functions / : K.^ ^ M 
which are continuous with period £ in each variable), via the maximum principle 
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(see e.g., [24]). Indeed, since is continuously embedded into Cj)(M^), the re- 
striction (^~^)[(7|j(ij2) of the operator ^q~^ to C)j(R^) is compact from Cj((M^) into 
itself. Moreover, it is clear that ^q"^ and (j</q^^)|c'j(i{2') have the same eigenvalues. 
Let now / be a non-negative (non trivial) function in Cjj(R^). Then, the function 
is in H^. Hence, in particular, it belongs to H'^{{—£,2i) x {—£,2i)) and 
solves the equation {fi + Xo)u — dAu = /. By the classical maximum principle, u 
is non-negative in {—£,2£) x {—£,2£). Actually u is everywhere positive. Indeed 
if u{xo) = at some point xq G {—£,21) x (— £, 2£), then, still by the maximum 
principle, it would follow that u = in [—£,21] x [—l,2£], which clearly cannot be 
the case. □ 

Appendix B. A maximum principle. 

Proposition B.l. Let be a second order operator with constant coefficients. 
Let u G satisfy the inequality J£u < 0. Then, u > 0. Similarly, if u belongs to 
C\{0,T);C{Tre)) nC{[0,T];H^), is such that ifu e C([0,T] x Qe), and it satisfies 
the differential inequalities Dtu — J/fu > and u < in [0,T] x dCli, then, u < 
m [0,T] X ile. 

Proof. For the reader's convenience, we sketch the proof of the second statement, 
the first one being a particular case of the second one. Since continuously 
embeds in the set of all continuous functions which are periodic, with period £ with 
respect to all the variables, then u can be extended by periodicity with a function 
(still denoted by u) which is continuous in [0, T] x and is therein continuously 
diffcrentiable with respect to the time variable. 

Suppose by contradiction that u is not everywhere non-positive in [0,r] x M.^ . 
Then, u has a negative minimum at some point {to,xo) € (0,To] x {-£,2£)^. Then, 
clearly, Dtu{tQ, xq) < 0. Moreover, since Jifu is a continuous function, ^u(tQ, xq) > 
0. The classical maximum principle yields the assertion. □ 
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